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Abstract 

Passivating lithium ion battery electrode surfaces to prevent electrolyte decomposition is crit- 
ical for battery operations. Recent work on conformal atomic layer deposition (ALD) coating 
of anodes and cathodes has shown significant technological promise. ALD further provides well- 
characterized model platforms for understanding electrolyte decomposition initiated by electron 
tunneling through a passivating layer. First principles calculations reveal two regimes of electron 
transfer to adsorbed ethylene carbonate molecules (EC, a main component of commercial elec- 
trolyte) depending on whether the electrode is alumina-coated. On bare Li metal electrode sur- 
faces, EC accepts electrons and decomposes within picoseconds. In contrast, constrained density 
functional theory calculations in an ultra-high vacuum setting show that, with the oxide coating, e~ 
tunneling to the adsorbed EC falls within the non-adiabatic regime. Here the molecular reorgani- 
zation energy, computed in the harmonic approximation, plays a key role in slowing down electron 
transfer. Ah initio molecular dynamics simulations conducted at liquid EC-electrode interfaces are 
consistent with the view that reactions and electron transfer occur right at the interface. Micro- 
gravimetric measurements demonstrate that the ALD coating decreases electrolyte decomposition 
and corroborate the theoretical predictions. 
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I. INTRODUCTION 



Improving the fundamental scientific understanding of lithium ion batteries^ - - is criti- 
cal for electric vehicles and other energy storage technologies. A key feature that enables 
the use of negative electrodes (graphite, Li metal, Si, Sn) operating below the reduction 
voltage of current commercial electrolytes is the formation of an electronically passivating 
but Li + -conducting solid electrolyte interphase (SEI) film on electrode surfaces.- - - Battery 
performance, irreversible capacity "loss," power fade, durability, exfoliation of graphite, and 
safety are highly dependent on the quality of the SEI. Therefore understanding the nature, 
formation composition, structure, and property of SEI is of great interest for Li-ion bat- 
teries. In this work, we apply computational and experimental techniques to analyze the 
success of the conformal atomic layer deposition (ALD) strategy for creating a passivating 
layer ( "artificial SEI" ) on electrodes,- - — focusing on graphitic carbon anodes.— 

It is generally accepted that, upon the first charge of uncoated graphitic anodes, the 
negative potential applied to induce Li + intercalation into graphite decomposes ethylene 
carbonate (EC) molecules in the solvent, yielding a self-limiting, 3-10 nm thick, passivat- 
ing SEI layer containing Li 2 C0 3 , lithium ethylene dicarbonate ((CH2C0 3 Li) 2 ),-^ and salt 
decomposition products. Early modeling work on organic solvent breakdown has focused 
on reactions inside bulk liquid regions, with an excess electron already injected.— - — While 
providing extremely useful predictions pertinent to that regime, such models necessarily 
ignore the possibility of surface-assisted reactions and effects arising from electron trans- 
fer from electrodes. A more rigorous if costly technique, ab initio molecular dynamics 
(AIMD), has recently been applied to simulate chemical reactions at several explicit solid- 
liquid interfaces.— - — One of the authors' previous AIMD works follows chemical reactions 
in real time at the pristine graphitic anodes/liquid EC interface.— 1 ^ It is found that, at the 
initial stage of SEI formation, fast e~ transfer and kinetically-controlled EC electrochemical 
reactions occur to form either CO or C2H4 gas,— - — mostly right at the oxidized edges of 
graphite sheets.— 1 ^ 

As electrolyte decomposition proceeds, e~ transfer becomes impeded by the intervening 
and partially-formed SEI layer between the solvent and electrode, and the decomposed 
solvent fragments can no longer anchor directly to the pristine electrode surface. This 
important next stage should figure equally prominently in the overall solvent breakdown 
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mechanism and the structure of naturally-formed SEI. The electron tunneling blockage by 
SEI layers is a kinetic (not thermodynamic) phenomenon, akin to stoppage of electron 
leakage through gate-oxide dielectric in semiconductor devices.— Theoretical study there is 
hindered by the substantial thickness, possibly porous/gel-like nature, and heterogeneous 
composition of natural SEI.— ~— i22 

Recently, it has been demonstrated that a sub-nanometer thick alumina layer created by 
the conformal ALD technique on graphite drastically diminishes solvent decomposition but 
permits lithium ion transport.— This promising ALD strategy even enables the cycling of 
low- melting-point propylene carbonate (PC), which otherwise exfoliates and destroys un- 
coated graphitic anodes. The mechanism of this ALD electrode passivation has not been 
completely understood. While expected to block or slow down electron transfer from the elec- 
trode to the solvent, it also appears to enhance the mechanical properties of the electrodes,- 
and likely hinders solvent intercalation between graphite sheets, thus preventing exfoliation. 
Apart from the technological implications, the unprecedented control over coating thickness 
and chemistry means that the ALD strategy also provides robust platforms for basic science 
studies of interfacial solvent decomposition reaction mechanisms, and for electron tunneling 
through the insulating layer which is a pre-requisite for electrolyte breakdown. 

In particular, the extreme thinness of ALD coatings lends itself to the present, predomi- 
nantly first principles computational study of electrode/solvent interfaces comprising up to 
850 atoms. Using crystalline, hydroxylated LiA102 layers as models of ALD coatings, we 
apply DFT and related theoretical techniques to show that ALD oxide films yield varying 
energetic, kinetic, and electron-tunneling impedence towards EC breakdown depending on 
the surface Li content and oxide thickness. We also provide evidence that e~ transfer occurs 
to EC molecules immediately next to electrode surfaces. Because these molecules are deep 
within the electric double layer (EDL), screening of electric fields by the EDL is less effective. 
In this sense, EC decomposition at battery anodes can differ fundamentally from classic elec- 
trochemical redox paradigms, where well-solvated transition metal complexes are separated 
by several Angstroms from the electrodes and "outer-shell" e~ tunneling dominates.— The 
solvent decomposition processes on ALD coatings provide insights that may be extrapolated 
to other passivation strategies, including natural SEI formed from electrolyte breakdown. 

Two limiting regimes of electron transfer, and two corresponding computational meth- 
ods, are emphasized. Rigorously, DFT deals with the electronic ground state, with nuclear 
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trajectories "adiabatic" to electronic configurations (i.e., ionic motions are slow compared 
with electron transfer). EC breakdown on Li (100) metal surface is in this adiabatic regime. 
DFT should be adequate for such processes, provided that the desired electronic configura- 
tion is the ground state and the self-interaction error of the approximate functional used is 
not critical to the properties being investigated.—"— 

In the opposite, non-adiabatic regime,—"— e~ transfer or tunneling is slow on the time 
scale of nuclear motion, and one must keep track of two electronic surfaces.— The oxide- 
coated model electrodes considered in this work pertains to this latter limit, where the 
electron transfer rate between two discrete orbitals is given by^ 



x/^ab| 
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A is the reorganization (free) energy, Vab is the coupling matrix element connecting the 
two electronic surfaces, and AG Q is the reaction free energy. A indicates the energy cost 
associated with molecular deformation needed to take on an extra electron (EC — > EC - ). 
Vab is the familiar prefactor that depends on the overlap between two many-body wave- 
functions associated with the two electronic surfaces (Fig. UK). Small Vab correlates with 
non-adiabatic e~ tunneling. 

Neither Vab nor A can be directly obtained using standard DFT methods. In this work, A 
is estimated using the constrained DFT (cDFT) approach^ 1 ^ and Marcus theory harmonic 
construction (Fig. [TJi) under both ultra-high vacuum (UHV) and liquid state^ 1 ^ configura- 
tions. cDFT is also applied to estimate Vab-— • While cDFT and related methods have been 
applied to molecules on metal surfaces,— — calculating Vab between a metallic electrode 
and an e~-accepting molecule, or for that matter the total e~ tunneling rate, has relied on 
simplified models.—"— When augmented using a Fermi Golden-rule expression (Fig. [lb), we 
argue that the our Vab value yields a well-defined kinetic prefactor for electron transfer from 
a metallic electrode. Our prefactor prediction is a preliminary estimate, and fundamental 
studies to extend cDFT to e~ transfer from metallic electrodes are needed. However, this 
is sufficient for our goal of order-of-magnitude estimates of e~ transfer rates. When the 
insulating layer (ALD oxide or natural SEI, or their combination) grows thicker, Vab starts 
to decay with oxide thickness, and its magnitude is examined via extrapolation in a way 
analogous to the one- dimensional Wentzel-Kramers-Brillouin (WKB) formula. An alterna- 
tive to this cDFT formulation may be Greens function/time-dependent DFT.—"— To our 



knowledge, TDDFT methods have not been successfully applied to predict orbital-to-orbital 
Vab values that involve metallic electrodes. 

With these computational techniques, we show that the sub-nanometer oxide coating,-^ 
generally not considered sufficiently thick for complete electron blockage in, say, gate oxide 
dielectric applications,— causes A (much neglected in previous battery studies) to play a 
significant role in ALD-assisted passivation. Electron tunneling to EC, not bond-breaking 
within the adsorbed molecule, is generally found to be the rate-determining step for break- 
down of EC adsorbed on the ALD-coated electrode. 

In terms of experiments, microgravimetric measurements that confirm the presence of 
solvent decomposition products on the surface are presented to corroborate aspects of our 
predictions. 

This paper is organized as follows. Section [XT] describes the methods used. Adiabatic 
electron-transfer induced EC reactions with Li metal surfaces is discussed in Sec. IHH The 
long-range electron transfer formalism is shown to be inapplicable here. Section [IV] de- 
scribes the non-adiabatic electron tunneling from oxide-coated electrodes to EC molecules 
adsorbed on their surfaces, and addresses the subsequent EC bond-breaking events. Adia- 
batic DFT/PBE calculations are shown to underestimate the electron tunneling barrier in 
this regime. Section |V] reports the experimental results, and Sec. |VT]briefly summarizes this 
work. 

II. METHODS 

A. Model systems 

The casual reader is encouraged to skip forward to Sec. [Ill] for the results. 

The key model systems are ~7.0 A ("thin") and ~10 A ("thick") layers of LiA102 
in /3-NaA102 structure with (100) surface terminations, coated on narrow strips of Li x Cq 
electrodes (Table [fl Fig. |2] more details of this oxide phase is provided in the supporting 
information (S.I.) document). Undercoordinated Al atoms on outer surfaces are terminated 
with OH groups, ensuring that surface states are removed. The oxide thickness is measured 
from Al to Al and excludes the surface hydroxyl groups or the C=0 edge atoms originally re- 
siding on Li x Cg. Crystalline LiA102 is a solid-state electrolyte candidate material.— LiA102 
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is used instead of AI2O3 to cover the possibility that the native AI2O3 layers deposited dur- 
ing ALD may have incorporated Li ions during the first charging half cycle. For example, 
some AlOH groups may be deprotonated at low voltages, causing Li + to coordinate to the 
A10~ and become part of the surface. The LiA102 mixed oxide thus allows us to exam- 
ine surface composition effects on EC breakdown. The stochiometry of the coating is such 
that their formal charges sum to zero. Another research group has found LiA102 signature 
on the surface of 5 nm thick ALD oxide films on Si anodes after power cycling using X- 
ray photo-electron spectroscopy. - Further computational evidence for Li-incorporation into 
AI2O3 films is presented in the S.I. If such Li + incorporation indeed occurs, the ALD layer 
will expand beyond its original AI2O3 thickness. 

The simulation cell, which provides a modest system size for AIMD simulations, is chosen 
so that the oxide is fairly well-matched to the Li x Cq surface cell, with oxide compressive 
strains of 1.8 % and 5.7% in the two lateral dimensions. The crystalline models are idealized; 
as in gate-oxide dielectric materials, insulating oxides should be amorphous to minimize 
cracks. The amount of Li present in the graphite region is determined by tuning the Li 
chemical potential to 2.1 eV. Upon geometry optimization, Li ions initially residing at the 
C=0 edges become strongly coordinated to the bottom surface of the oxide coating. Another 
model, with a single 10 A thick layer of LiA102 hydroxylated on both sides but no Li^Cg 
component, is used to examine post e~-transfer EC - bond-breaking. 

To emphasize the influence of surface groups, we also include a model with a ~5.0 A 
thick layer of a-A^O^ coated on both sides of the Li x Cq strip. The oxide layers have (0001) 
terminations with AlOH surface groups (Table [I]). It has been predicted that 7-AI2O3 is 
more stable than the a phase for film thickness below 36 A.— However, this estimate was 
made without accounting for surface hydroxylation. Since our thin AI2O3 film contains only 
two Al-0 layers (not counting the C=0 edge groups), the oxygen positions are arguably 
consistent with both a-Al 2 03 with closed packed oxygen in ABAB stacking, and cubic 7- 
AI2O3 with ABCABC stacking. In a-Al 2 03, all Al are in octahedral sites while Al occupy 
both octohedral and tetrahedral sites in 7-AI2O3. Upon applying geometry optimization to 
the initial "a-ALjCV film, some Al ions are found to migrate to tetrahedral sites, especially 
those coordinated to graphite-edge C=0 groups. Thus our "a-Al 2 3 " film arguably exhibits 
both a and 7 character, consistent with experimental ALD coatings which are considered 
amorphous without long-range order. 
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system/coating 


method 


Figure 


stochiometry 


cell size 


N EC 


thin LiA10 2 


NEB 


Fig. 


Hb-d,E 


Al4 8 Li 96 Oi4§ Cg 2 H 24 


40.00x12.47x15.06 


1 


thick LiA10 2 


NEB 


NA 




Al7 2 Liii8O 2 08Cg 2 H 2 4 


46.0x12.47x15.06 


1 


thick LiA10 2 


AIMD 


Fig. 


m 


Al7 2 LingO 2 08Cg 2 H 2 4 


48.5x12.47x15.06 


36 


thin LiA10 2 


AIMD 


Fig. 




Al48Lig6 Oi48 Cg 2 H 2 4 


43.00x12.47x15.06 


36 


only LiA10 2 


NEB 


NA 


Al36Li360s4H 2 4 


24.00x12.47x15.06 


1 


A1 2 3 


AIMD 


Fig. 


m 


Al7 2 2 o4Ci 2 oH7 2 Li5i 


33.34x14.97x18.82 


36 


Li (100) 


AIMD 


Fig. 




Lig6 


30.35x14.63x14.63 


32 


Li (100) 


NEB 


Fig. 




Lig6 


30.35x14.63x14.63 


1 


Li (100) 


NEB 


Fig. 


[3b-d 


Li48 


24.00x9.75x9.75 


1 



TABLE I: Details of systems used in AIMD and geometry optimization-plus-NEB barrier calcu- 
lations. The spatial dimensions are in A. "Stochiometry" omits ethylene carbonate atoms in the 
liquid region. 

Finally, a thin slab of lithium metal truncated along (100) surfaces is considered. Even 
though Li metal itself cannot currently be used as rechargeable anodes, EC breakdown 
products on Li are qualitatively similar to those on LiC6 surfaces.- 1 ^ Under open circuit 
conditions, Li metal should be at a well denned ~ —3 V versus the standard hydrogen 
potential.— Furthermore, EC decomposition on Li surface is free of the ambiguity associated 
with solvent co-intercalation into graphite.— Thus Li metal provides an useful baseline with 
which to interpret predictions for the oxide-coated surfaces. 

B. Adiabatic regime: DFT, AIMD simulations 

All calculations are performed using the Vienna Atomic Simulation Package version 4.6 
(VASP)^ 1 ^ and the PBE functional.— AIMD simulations apply T-point Brillouin zone sam- 
pling, a 400 eV planewave energy cutoff, and a 10 -5 eV or 10~ 6 eV convergence criterion 
at each Born-Oppenheimer time step. The trajectories are kept at an average tempera- 
ture of T=450 K using Nose thermostats, except for the EC/Li metal simulation where 
T=350 K is enforced. Tritium masses are substituted for protons to enable a time step of 
1 fs. Under these conditions, the trajectories exhibit drifts of less than 1 K/ps. Due to the 
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approximate nature of DFT functionals and the simulation protocol (tritium masses and 
thermostat used), the predicted reaction time scales should be treated as relative, not abso- 
lute. AIMD simulations reported do not account for spin-polarization. Our previous work 
has revealed no qualitative difference between restricted singlet and spin-triplet DFT/ AIMD 
simulations.— Molecular configurations are pre-equilibrated using Monte Carlo simulations 
and simple molecular force fields, as described in an earlier work.— Representative AIMD 
snapshots are depicted in Fig. |2j 

The AIMD liquid/solid interfacial simulations are akin to dipping electrodes fully pre- 
intercalated with Li into the organic solvent. In principle, it may be possible to intercalate 
Li + in the electrolyte, remove the anodes from solution, clean off possible decomposition 
products in inert environments, and re-insert in solution to measure the open circuit voltage. 
Such experiments have not been performed but can be attempted in the future. In Sec. IIVDI 
we further discuss the electrochemical potential of these electrode models. 

T=0 K geometry optimizations and climbing image nudged elastic band (NEB)^ bar- 
rier calculations (e.g., Fig. [3]) are performed with spin-polarization, a 10 -4 eV convergence 
criterion, and a linear potential correction applied in the direction perpendicular to the sur- 
face to remove dipole-image interactions.— T-point sampling is generally applied, except for 
calculations involving Li metal slabs where 1x2x2 Brillouin zone sampling is used. Even 
there, T-point NEB calculations yield a C-0 bond-breaking barrier only 0.1 eV higher than 
the more dense Brillouin grid result. It is also found that the geometry and net charge of an 
adsorbed, intact EC - on LiA102 surface is unchanged whether T-point or 1x2x2 grids are 
used. Comparing restricted singlet and spin-polarization results, no difference is discernable 
in the EC on Li metal calculations, where the bond-breaking barrier is small (< 0.1 eV, 
Sec. IIIII) and adiabatic electron transfer from the electrode and the bond-breaking event 
occur simultaneously. These are the conditions under which spontaneously EC decomposi- 
tions are observed in picosecond AIMD simulations, justifying the use of non-spin-polarized 
DFT there. Higher bond-breaking barriers, like those on the 10 A thick LiA102 surface 
(see the S.I.), are reduced when spin polarization is allowed. A spot check shows that spin- 
unrestricted DFT calculations reduce the Cc-0 cleavage barrier by 0.15 eV on this surface. 
Even with this reduction, the barrier is high enough to prevent observation of EC breakdown 
in picosecond time scale, and therefore using non-spin-polarized DFT in AIMD simulations 
does not affect the conclusion that no reactions occur within the 7 ps trajectories in high 
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barrier cases. Further details on NEB calculations are discussed in the S.I. 
C. Non-adiabatic regime: Constrained DFT 

A version of the constrained DFT (cDFT) method^^ is implemented into VASP. The 
constraining potential is chosen to be 

W(v) = Voll-IkMr)], (2) 
fi(r) = [1 + tanh(«(|r - r<| - Wi )]/2, (3) 

Here V is a constant to be self-consistently determined, k is 6A _1 , % labels the atoms in the 
selected EC participating in electron transfer, is the atom position on that EC, and Wi is an 
element-specific radius. Wi amounts to 1.65 A for C and O and 1.25 A for H. These values are 
similar to Lennard- Jones radii in simple atomic force fields. A more stringent wavefunction 
convergence criterion of 10 -6 eV or smaller is enforced in self-consistent cDFT calculations. 
The W(r) functional form does not double-count electron density on adjacent atoms and 
appears pertinent when bond-breaking can occur. Normalized, atomic orbital-based charge 
projection operators used in the literature^^ 1 ^ may be less applicable for electron transfer 
coupled to bond-breaking, but they can be tested for the present application in the future. 

The total electronic charge on the selected EC is determined by projecting W(r)/V on to 
the DFT electron density. With the w% values mentioned above, unconstrained DFT predicts 
that a charge-neutral EC molecule adsorbed on the thin LIAIO2 surface (Fig. Hb) exhibits a 
slight +0.20 1 e | "net charge," while — 0.60±0.1|e| resides on the EC~ (Fig. HH). The adsorbed 
EC - exhibits a similar — 0.67|e| charge on the L1A102 oxide slab without any conductive 
Li x C6 component (Tabled]). The non- integer values arise because of residual charge densities 
at the edge of EC molecules beyond the range of W(r)/V . (The net spin on EC~ is about 
0.9|e|, and is more centered on EC than the net charge.) Increasing is ruled out because 
of the close proximity of adsorbed EC to the surface hydroxyl groups. For example, using 
larger Wi has been found to lead to abstraction of protons from surface hydroxyls. The 
protons then bind to the negatively charged EC molecule. Such reactions are not seen in 
unconstrained AIMD simulations and are deemed unphysical. We have therefore defined 
+0.20|e| and — 0.60|e| to be the net charges of fiat, intact EC (Fig. Hb) and EC" (Fig. HH) 
when using self-consistent cDFT calculations to impose charges on the molecule. Increasing 
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\V \ to increase the charge on adsorbed EC" to — 0.80 |e| is found to yield only a 10% change 
on the coupling matrix element Vab, but can increase A by a fraction of an electron volt. The 
more important parameter, the barrier in Eq. (TJ is only affected by ~ 5X/A in the harmonic 
approximation used in this work. In the S.I., the predicted A for adsorbed EC is shown to 
be comparable to that for EC in liquid EC, computed using cluster calculations, localized 
orbitals, and a dielectric continuum approximation. 

Coupling matrix elements Vab between the two different adiabatic surfaces (Fig. [TJi) are 
computed using the cDFT formalism for discrete orbital levels,— which is implemented into 
VASP within the projector-augmented wave formalism.— The same atomic configuration 
must be used for both electronic surfaces, and this is chosen to be the optimized atomic con- 
figuration where no excess electron resides on the flat, adsorbed EC. Vab is generally assumed 
to be relatively independent of atomic positions with the "Franck- Condon" approximation, 
although molecular orientation dependence has been demonstrated.— Vab emerges from the 
2x2 Hamiltonian matrix H connecting the donor i n our case from unconstrained 

DFT calculations) and acceptor (|$b)) single determinantal wavef unctions.— 1 ^ |$b) fea- 
tures an excess electron on one EC molecule and is generated using cDFT. H contains the 
overlap matrix element as well as ($a\ S e W(r e )\§B) , where e labels all occupied 

electronic levels.— These calculations are fairly costly and are performed at T=0 K in this 
work. 

This cDFT-based Vab formulation was originally devised for electron transfer between 
ground state cDFT donar and acceptor electronic configurations, with the implicit assump- 
tion that the relevant density-of-state is discrete. In the limit of non-interacting electrons 
residing on a metal electrode, this formalism reflects only the top curve on the left side of 
Fig. [lb and does not reduce to the well-known Fermi Golden Rule formula for tunneling 
from a continuum of donor states. Consider, in this limit, a band of single-particle energy 
levels E characterized by a density-of-state D(E) of orbitals 4>{E), Fermi distribution func- 
tion f(E), Fermi level Ep, and an isolated acceptor orbital <p a with energy E a . The Golden 
Rule rate, associated with multiple level crossings illustrated in Fig. [lb, is 

fc GR oc^y dE\( ( p(E)\v(E)\ ( j )a )\ 2 D(E)f(E-E F )5(E-E a )^ , (4) 

where v(E) is the single-particle coupling matrix element and (O)r denotes averaging over 
nuclear degrees of freedom R on which all quantities implicitly depend. This formula allows 
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many-electron acceptor \<$>b') states that involve a but not the HOMO of which rep- 
resent electron-hole excitations. 52 In contrast, cDFT can only generate the electron-acceptor 
manifold \<&b') which is the ground electronic states within the applied constraint. 

To incorporate the effect of Eq. HI we make the common assumption that Vab is con- 
stant over the relevant range of density-of- state.— >^~— Then an empirical Golden Rule-like 
expression can be proposed: 



/r GR _ y f V^I^ABl 2 



(AE Q + AE a , + A) 
4Xk B T 



(5) 



Here AE is used in place of AG Q because we ignore entropy changes in T=0 K, UHV-setting 
calculations, f a > is the Fermi and/or symmetry weight of Kohn-Sham orbital a', and AE a i 
is the difference in energy between the Fermi energy and each Kohn Sham orbital level a', 
e F —e a >. a' deep within the occupied manifold does not contribute due to the AE a > factor. 
The self-consistent T-point electronic density is used to generate a dense grid of occupied 
states <fi a i using a 1x4x4 Monkhorst-Pack Brillouin sampling.— 

To converge to the infinite size limit for e~ transfer to a single EC molecule, the correct 
approach is not to increase fc-point sampling, but to increase all spatial dimensions of the 
model electrode. If the Li x C6 component of the electrode is doubled in size in any one 
direction, the orbital donor wavefunction a / delocalized over the electrode is scaled down by 
~ I/a/2, and |Vab| 2 decreases 2- fold. This underscores the fact that Vab is not a measurable 
quantity in finite-sized electrode models, but changes with the system size. However, the 
density-of-state D(E a r) increases proportionately with system size, and the sum over all 
orbital contributions (Eq. [5]) should be well-defined in that infinite size limit. 



D. Experimental Details 

Carbon films deposited onto Cu were used as electrodes to explore the passivating role 
of the ALD-derived alumina coatings with respect to electrolyte reductive decomposition. 
Polished AT-cut quartz crystals patterned with Cu electrodes (9 MHz, Inficion) were used 
as the base current collector for conducting both voltammetry and gravimetry. 50 nm thick 
carbon films were deposited onto these crystals using a pulsed laser deposition method.— 1 ^ 
Conformal alumina coatings were deposited onto both carbon films and bare Cu electrodes 
at a substrate temperature of 180 °C using alternate cycles of trimethylaluminum and water 
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to produce amorphous AI2O3 films of either 0.55 or 1.1 nm thickness.— — iiZZ Cycles of NO2 
and TMA pre-exposure were used to ensure the nucleation and growth of a continuous 
alumina film. 7 - Electrochemical measurements were conducted under argon in a glove box 
(Vacuum Atmospheres, <100 ppb H 2 0, <1 ppm O2) in 1 M LiPF6 in a 1:1 volume mixture of 
ethylene carbonate and diethylcarbonate (Hoshimoto and Kishida Chemical). A Solartron 
1287 potentiostat coupled with a Maxtek RQCM controller were used for simultaneous 
voltammetric and gravimetric measurements. 



III. RESULTS: ADIABATIC AIMD/DFT PREDICTIONS OF EC/LI(100) REAC- 
TIONS 

Adiabatic DFT/PBE calculations should be pertinent for the EC/Li(100) interface, where 
EC and the metallic electrode are in close contact and fast e~ transfer is expected. 



A. Liquid EC on Li (100) 

Liquid EC has been previously predicted to decompose at the C=0 edges of LiCg elec- 
trodes within 7 ps at T=450 K in AIMD/PBE simulations.— This timescale is used to 
qualitatively gauge the DFT/PBE predicted reactivity of other surfaces towards liquid EC. 

Figure [2b shows that liquid EC decomposes readily on Li (100). Within 15 ps, all 12 
EC molecules adjacent to the Li metal, out of 32 EC in the simulation cell, have accepted 
electrons and decomposed. 11 out of these 12 exhibit two broken Cc-0 bonds to form CO + 
OC2H 4 2- ;— "— only one EC decomposes in the classic C2H4 + CO 2 " route hitherto widely 
accepted in the literature, cleaving both Ce-0 bonds. 1 ' 2 Here Cc and Ce are the carbonyl 
and ethylene carbon atoms, respectively. This finding is consistent with those in Ref. |21 



where both CO and CO 2 products emerge at the interface between liquid EC and pristine 
LiCe with oxidized edge groups. This agreement is significant because, by construction, 



the models used in Ref. 



21 



exclude solvent co-intercalation cited in the "3-dimensional" 
SEI formation pathway.™ Nevertheless, fast EC decomposition and identical products are 
predicted on both pristine graphite and Li metal surfaces, showing that such co-intercalation 
is not necessary for SEI initiation. 

In the EC/Li trajectory, the temperature is thermostat at T=350 K, not T=450 K, 
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the 



15 



to avoid melting the solid Li. Despite this, the heat generated by the reactions and 
incorporation of CO into the metal slab have caused significant amorphization. In Ref. 
the initial 200 fs of this trajectory is examined in detail. The bent EC geometry, with 
the carbonyl C=0 displacing out of the EC plane, is shown to be correlated with electron 
transfer to EC, just like for the isolated EC - in solution (Fig. [TH).— ^ This bent geometry 
plays a critical role in electron transfer and reorganization energy calculations in L1A102- 
coated surfaces (see below). Our AIMD simulations have shown that OC2H4O 2- can react 
with 2 CO2 to form the main SEI organic product ethylene dicarbonate. Whether this 
product is deposited at the initial stage of SEI growth depends on the availability of CO2 
and the solubility of the decomposition fragments-^ 



B. Isolated EC on Li (100) 

Remarkably, even a single EC molecule, in the absence of the liquid environment which 
stabilizes its ionic breakdown products, still decomposes on Li (100) surfaces to form CO + 
OC 2 H 4 2 ~ within picoseconds. (Fig. [2b ) This suggests that a simple T=0 K energy profile 
calculation is relevant to EC decomposition.— 

Figure [3b, compares the T=0 K energy profiles of the two modes of excess electron- 
induced EC breakdown on Li metal. They show that cleaving the Cc-0 bond to form the 
precursor to carbon monoxide, OC2H 4 OC0 2 ~, is thermodynamically less favorable than the 
ethylene carbon-oxygen bonds to form CO 2- and C2H4 by a substantial 1.53 eV. Cleaving 
the remaining Cc-0 bond in the CO route only leads to another 0.16 eV stabilization. The 
barriers associated with both types of bond-breaking are vanishingly small. Applying the 
HSE06 truncated hybrid functional,—^ 2 - which exhibits far less self-interaction errors^ 2 - - — 
than PBE, increases the Cc-0 breaking barrier, but only to 0.16 eV (not shown, but consis- 
tent with the similar short time dynamics predicted with the PBE and HSE06 functional^) . 
This suggests that adiabatic DFT/PBE barrier predictions are reasonably accurate for EC 
in contact with Li metal. The small barrier explains why both product channels are available 



in picosecond time scales at explicit liquid EC/electrode interfaces (Fig. |2b, Ref. Il5|). We 
speculate that the kinetic prefactor favors the CO-route and makes it the majority product 
in liquid-solid interface simulations (Fig. |2b). 
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C. Long-range e transfer formalism is not applicable to EC/Li(100) 

For e~ transfer to EC directly adsorbed on uncoated electrode surfaces, the close contact 
should render the cDFT method for non-adiabatic long-range electron transfer— ^2 inappli- 
cable. If one insists on calculating Vab using cDFT and and the simulation cell described in 
Table [U Vab is found to be 0.23 eV for a flat EC adsorbed on Li (Fig. [3b). This large Vab is 
consistent with the significant, 56% overlap between the acceptor and donor many-electron 
wavefunctions, and should put the system in the adiabatic electron transfer regime — even 
with the caveat about the system size dependence of Vab-— (For comparison, a theoretical 
work on NO molecules adsorbed on Ag(lll), not using cDFT, has also yielded fraction-of-eV 
Vab-—) We conclude that the adiabatic DFT/PBE treatment should suffice in this case. 

IV. RESULTS: NON-ADIABATIC ELECTRON TRANSFER TO EC ON OXIDE 
SURFACES 

This section focuses on a UHV-like model consisting of an isolated EC adsorbed on the 
lithium- intercalated graphitic carbon strip coated with LiA102- A 0.4 V/A electric field 
is applied. For this model, e~ tunneling resides in the non-adiabatic regime where cDFT 
calculations are pertinent. The relevance of this model to the liquid EC/electrode interfacial 
environment will be clarified below. 

A. Two metastable EC charge states on 7 A thick oxide surface 

The Lia-Cg model with a 7 A thick LiA102 coating proves especially useful for examining 
the details of electron transfer from the electrode to an adsorbed EC, which either precedes 
or takes place simultaneously with EC~ decomposition. Two (meta)-stable adsorbed EC 
configurations can be stabilized (Fig. Hj). One is a flat, charge-neutral EC coordinated to a 
surface site (an AlOH group) via its carbonyl oxygen atom (Fig. Hb). Figure depicts the 
local electronic density-of-state (DOS) for this system. The Li x C6 region contains partially 
occupied states near the Fermi level (E F ). The insulating oxide spans a substantial band 
gap, although there are surface states in the interface with Li^Cg that reduce the effective in- 
sulating thickness. The highest occupied molecular orbital (HOMO) of EC is below —2.5 eV 
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while the LUMO lies above E-p. This DOS is consistent with a charge-neutral EC weakly 
interacting with the oxide surface. 

The other configuration has an intact EC - which adopts a bent geometry with the C=0 
bond protruding out of the EC plane (Fig. HU). This is reminiscent of the first stage of 
liquid EC decomposition on Li (100) surface, where the e~-accepting EC adopts a similar 
bent configuration.— The excess charge on the EC is centered around the carbonyl oxygen 
atom which is coordinated to two AlOH groups and a Li surface atom. The system exhibits a 
DOS (Fig. Hb) substantially different from Fig.0|i. The majority spin, highest occupied state 
of the EC molecule now lies below the Fermi level. The shift in the LUMO upon e~ addition 
serves as a caveat against using the LUMO of the neutral molecule as a figure-of-merit in 
assessing electrochemical reduction tendencies. 

The bent EC~ is almost iso-energetic with the flat EC. Its slight exothermicity, AE Q = 
—0.02 eV, does not depend on whether the electron transfer is adiabatic or non-adiabatic. 
It should not be affected by the periodic images imposed by the simulation cell because the 
dipole correction is applied.— In fact, despite the transfer of an e~ across a 7- or 10-A thick 
oxide layer, the overall dipole moment of the simulation cell changes by less than 1.0 |e|A, 
apparently because the electron density in the metallic Li^Ce strip can rearrange itself to 
accommodate the electron transfer. The total charge in the simulation cell is conserved in 
these calculations and the large correction due to periodic boundary conditions for isolated 
ions in solutions is not needed.— 1 ^ Note that AE Q is used in place of AG a because the 
calculation is performed at T=0 K. 

B. Non-adiabatic electron transfer on oxide surface 

We apply the cDFT method to calculate A and Vab required to estimate the electron 
transfer rate k et (Eqs. [T]&[5]). We stress that the flat EC absorbed on the oxide coatings is 
treated using unconstrained DFT/PBE. The highest-occupied orbitals of the 7 A and 10 A 
thick coatings reside in the Li x Cg region, and exhibit integrated electron densities of less 
than 10~ 4 and 5xl0 -8 |e| on the EC molecule, respectively. This shows that the uncon- 
strained DFT method already gives a reasonable description of the neutral EC electronic 
configuration. 

A is computed for the optimized, flat EC geometry adsorbed on the thin LiA102 coating 
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(i.e. image in Fig. [5^). cDFT imposes an extra electron on the EC molecule. On the 
7 A thick coating, it yields a vertical excitation energy AE vevt =X + Ai? o =2.04 eV, where 
AE Q is the aforementioned —0.02 eV offset between donor and acceptor. Alternatively, an 
electron can be removed from the frozen bent EC - configuration (image 5), which leads 
to A' + A£' o =1.80 eV. A and A' agree to within 14%. This is qualitatively consistent with 
the Marcus theory postulate that the polarization degrees of freedom respond harmonically 
(Eq. [1]), yielding a single reorganization energy that governs electron transfer reactions.— 
With A=2.06 eV for EC adsorbed on the thin LiA102 surface, the non-adiabatic barrier 
becomes 0.51 eV from a simple Marcus construction (Eq. [1]). This barrier is much higher 
than the ~0.1 eV adiabatic DFT/PBE activation energies for both the Cc-0 and Ce-0 
bond breaking pathways on this surface (Sec. IIVCI) . and is therefore the rate-limiting step 
in EC breakdown on the surface of the thin LIAIO2 coating. 

In the S.I., an EC with a dielectric approximation of the liquid EC solvent medium is 
found to exhibit an average of A = 1.76 eV, similar to EC adsorbed on the thin LIAIO2 
coating. The co-solvent dimethyl carbonate (DMC) exhibits only slightly smaller A values. 
Therefore the substantial A, large compared to many organic molecules,— is intrinsic to 
out-of-plane bending of the C=0 group as the carbonyl carbon atom adopts a sp 3 -like 
geometry to accommodate an e~. The S.I. further presents results on vertical excitation 
energy, A^V^t = A + AE D (Fig. UK), computed in several AIMD snapshots, to suggest that 
the Arrhenius term in Eq. [JJ favors e~ transfer to EC molecules at the interface over EC in 
the bulk liquid region. In such AIMD simulations, we are limited to the first choice of A, 
i.e., instantaneously adding an electron to EC, because EC - in liquid EC can have short 
lifetimes.™ Hence we will focus on this first choice throughout this work. 

As this is a T=0 K calculation in a UHV-like setting, we have simply used the (AE a + 
A) 2 /(4A) expression in Eq. [1] as the tunneling barrier,— ^2 and have not traced out the 
two adiabatic curves as a function of the energy gap using liquid state potential-of-mean- 
force simulations.— 1 ^ We have however checked that, when relaxing EC - frozen in the flat 
geometry (Fig. Hb) with a constrained charge, it reverts to the stable bent EC - (Fig. HU) 
configuration, showing that the cDFT approach puts the system on the correct electron- 
acceptor potential surface. In the future, we plan to perform direct cDFT calculation of 
the barrier height at T=0 K by simultaneously optimizing the same atomic configuration 
on both energy surfaces. 
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coated electrode 


uncoated electrode 


system 


work func. 


system 


work func. 


thin LiA10 2 (OH) 


2.47 


Li(100) 


3.05 


thick LiA10 2 (OH) 


2.90 


graphite edge 


4.57 


thin LiA10 2 (OLi) 


2.25 


LiA10 2 


5.42 


A1203 (OH) 


4.10 


A1 2 3 (0001) 


6.22 



TABLE II: Work function of model systems used in this work computed using the PBE functional, 
in eV. The left column describes the oxide coatings on LLjCe; the right column refers to work 
functions of pure crystals. The Al 2 03 (0001) model is Al-terminated. The graphite slab has 
dangling bonds. 

The cDFT coupling matrix element is estimated to be Vab=0.022 eV at the flat EC 
geometry. Fig. [5fc depicts the highest-occupied DFT and cDFT orbitals, integrated over 
the lateral dimensions, for the systems with flat EC and flat EC~ respectively. The overlap 
between them, 

(0ec M °I0ec-°)> is 0.0125, or within 2% of that between the respective 
determinantal wavefunctions (^aI^b) (Sec. [XT]) which includes many-electron contributions. 
Therefore the relaxation of other electrons ( "polarization effect" ) does not strongly influence 
the overlap integral when using T-point sampling. 

This estimate of Vab does not reflect the classic Fermi Golden rule phenomenology 
(Sec. Hi])- Applying Eq. [5] to approximately account for the finite density-of-state on the 
electrode, we obtain a 1.63xl0 4 /s electron transfer rate. Simply using the cDFT definition 
of Vab in Eq. [TJ which represents a single point integration quadrature, merely underesti- 
mates this rate by a factor of 1.68. Using DFT/PBE rather than more accurate but costly 
hybrid functionals has been known to overestimate Vab by almost a factor of 10.— In the 
present case, the DFT/PBE underestimation of the band gap of the insulating oxide layer 
may lead to some overestimation of the electron tunneling rate. Despite the approximations 
and assumptions involved, this is to our knowledge the first DFT-based estimate of the tun- 
neling rates from an electrode, through an oxide layer, to an adsorbed EC molecule. The 
value may potentially be compared with UHV measurements. After electron transfer, EC - 
decomposes, and the negatively charged EC fragments will most likely complex with Li + 
from the electrolyte and be incorporated into the SEI layer on top of the ALD film. 
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On the thicker LiA10 2 coating, A£ o =0.77 eV. A=1.98 eV is predicted in the flat EC 
geometry. Removing an e~ from the bent geometry yields A' = 1.69 eV. The asymmetry 
is 15%. We again adopt the first choice of A. Vab is estimated at 0.0128 eV, about half 
that of the 7 A thick LiA102 coating. 9 - As discussed in Sec. IIV} the thinner coating exhibits 
substantial surface relaxation which is absent in the 10 A layer, making a purely thickness- 
based comparison of Vab difficult. Fig. [5U depicts the donor and acceptor Kohn-Sham and 
cDFT orbitals. The overlap between them, (^ecI^ec-)) is about 0.004, a factor of 3 less 
than that across the 7 A thick coating. Including the contributions of Eq. [5J k et becomes 
extremely small (2.8xl0 _5 /s) due to the larger AE a . 

The overall k et is clearly very sensitive to AE a or AG . In UHV settings, AE depends 
on both the electric field and surface heterogeneity at atomic lengscales (see below). At 
electrode/liquid electrolyte interfaces, k et is a function of the applied voltage as well as the 
local EC reduction potential via AG Q (Eq. [I]), which may be a function of the distance from 
the electrode. Direct measurement of the reduction potential of an intact EC is unavailable 
because EC decomposition occurs faster than cyclic voltammetry time scales. If one adopts 
a theoretical AG Q =— 0.15 eV for e~ transfer to intact EC molecules in EC liquid in at 
Li(s)/Li + voltages (S.I.), the predicted initial electron transfer rates k et through the 7- 
A and 10-A thick oxide coatings (~ 1.7xl0 5 /s and 8.3xl0 4 /s), will permit electrolyte 
breakdown, even if we assume that these rates are overestimated by 100 times due to the 
use of the PBE functional discussed above. Indeed, our gravimetric measurements reveal 
electrolyte decomposition on the coated electrodes - - consistent with ready availability 
of electrons - - albeit in much less quantity than on uncoated electrodes (Sec. |V|). The 
electrolyte decomposition product then yields an additional insulating layer that prevents 
further electron tunneling. 

Our main point in this section is not to predict exact k et values, but to highlight the 
previously neglected role of the EC reorganization energy (A) on electrode coated with an 
insulating layer. An immediate implication is that different solvent molecules/salt compo- 
nents may exhibit different A and e~-transfer rates. 
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C. DFT/PBE treatment of electron transfer on oxide surface is inadequate 



We next demonstrate that adiabatic DFT/PBE calculations are inadequate when dealing 
with e~ tunneling through insulating oxide layers. 

The electron transfer barrier strongly depends on whether the e~ transfer is adiabatic 
or not, and on the accuracy of the DFT method used. Figured depicts a climbing-image 
NEB calculation with 4 images along the reaction coordinate linking the flat EC and the 
bent EC - to examine the DFT/PBE adiabatic energy landscape in the 0.4 V/A electric 
field. DFT/PBE predicts a 0.09 eV barrier associated with electron transfer through the 
thin LiA102 layer. 

This small 0.09 eV value gives the strongest indication that DFT/PBE grossly underesti- 
mates the e~ transfer barrier. In classical electron transfer paradigm (Fig. UK), the parabolic 
intersection which yields the non-adiabatic barrier in the exponential term in Eq. [1] is ex- 
pected to differ from an adiabatic prediction of barrier by Vab- Instead, the former is 0.51 eV 
and the latter is 0.09 eV (Fig[5h); their difference far exceeds Vab=0.022 eV before con- 
sidering system size dependence. The discrepancy is most likely due to the self-interaction 
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The widely 



error in the DFT/PBE functional,— 1 ^ a point already alluded to in Ref. 
used PBE functional, along with others, do not sufficiently penalize configurations where an 
electron occupies both the electrode and the EC molecule. Indeed, in image 2 of Fig. [5^, a 
fractional —0.2 |e| charge develops on the EC, which should be considered unphysical for a 
molecule separated from the electrode by at least 7 A. Hybrid DFT functionals exhibit less 
self-interaction errors than DFT/PBE, but are currently too costly for computing barriers 
in interfacial systems of this size.— 

The 10 A -thick oxide-coated electrode exhibits a monotonic DFT/PBE energy profile for 
electron transfer. There is no DFT/PBE adiabatic barrier between the flat EC and bent 
EC~ beyond the minimal 0.77 eV mandated by the endothermicity (Fig. \5jp), suggesting 
that the electron tunneling barrier is again severely underestimated. Using the conjugate 
gradient geometry minimizer in VASP, the bent EC - geometry on this surface is in fact on 
the verge of instability, about to lose electron density to the electrode and relax to the flat 
EC geometry. Therefore the depicted energy profile actually reflects an optimized geometry 
subject to a charge constrained via cDFT with a small V = —0.2 eV. 
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D. Work function and electrochemical potential 



The electron tunneling rate at electrolyte-electrode interfaces depends on the electro- 
chemical potential (<&) of the electrode. In the coated graphite model systems, $ is not 
precisely known. Directly calculating $ involves averaging the electrostatic potential differ- 
ence between the conductive (inner) region of the electrodes and a distant point in the bulk 
liquid beyond the thickness of the electric double layer,— and involves consideration of image 
charge and surface potential effects.—"^— These are beyond the time and length scales of 
current AIMD simulations. Fortunately, the EC/Li(100) interface mimics immersing freshly 
prepared Li metal into liquid EC, and reflects an unambigous open-circuit voltage below the 
threshold at which EC becomes electrochemically decomposed (+0.8 V vs. Li + /Li(s)). This 
is a major reason Li is considered in this work. 

If we consider the energy of an e~ in the bulk electrolyte to be a constant, independent 
of electrode surfaces, the energy for ejecting an electron from different electrodes into the 
bulk electrolyte will only be shifted by the work function^ (where an e~ goes into vacuum). 
Thus, we have computed the work functions of coated and uncoated electrode surface and 
some crystal planes of ALD coating materials (Table [TTJ) . The -OH and -OLi terminated 
LiA102 coating work functions are within 0.5 eV of the Li metal value, indicating that 
similar energies are required to remove an electron from these surfaces. The AI2O3 coated 
surface has a much higher work function (Table HI]), consistent with our observation that 
AI2O3 is a more insulating material than LiA102 (see below). 

Even though our DFT calculations show that placing these oxides in contact with Li 
metal surfaces leads to immediate Li metal oxidation, we use Li(100) as a reference because 
its voltage is similar to that of LiC^. Aligning the work functions of Li(100) and the oxide 
materials (Table HI1). it is clear that the valence and conduction bands of the ALD phase lies 
below and above the Fermi energy (Ep) of Li metal, respectively. Electron tunneling from 
the Li Ep to the conduction bands of Al-terminated AI2O3 (0001) and LiA102 (100) exhibit 
1.43 eV and 1.13 eV offsets (barriers, AE), respectively. 

According to the ID WKB formula, the tunneling prefactor is 



where m e is the electron mass. If we take a tunneling transmission probability of e 40 as the 
limit of vanishing electron tunneling, 3.7nm thick LiA102 and 3.2nm thick AI2O3 are required 




(6) 
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to stop total SEI growth using DFT/PBE predicted AE. The work function is only one 
contribution to $ and does not contain solvent orientation and electric double layer effects^ 
(which should be less important for our inner-shell redox reduction of solvent compared to 
the classical paradigm of electron transfer to well-solvated outer-shell ions). Nevertheless, 
it gives a simple guidance for comparing different insulating ALD coating materials. As e~ 
transfer slows down and becomes rate-limiting, the composition of SEI films formed from 
electrolyte decomposition will likely change. This is because solvent molecules (other than 
EC), the counter ions (PFg) in the salt, and other partially decomposed products may 
exhibit smaller electron transfer barriers (reorganization energies) and start dominating the 
product channel. 

E. EC bond-breaking on ALD coating after e transfer 

On the 10 A thick LiA102-coated Li x Cq strip (Fig. [2]F), no EC decomposes within 
7 ps. The limited duration of the AIMD trajectory does not permit an estimate of the 
adiabatic AIMD/PBE free energy barrier. While this barrier can be computed using 
the AIMD/potential-of- mean- force method-^ 7 it will be underestimated due to PBE self- 
interaction errors and underestimation of the electron tunneling barrier. 

However, on the 7 A thick LiA1022 layer, a Cc-0 bond on one EC molecule is sponta- 
neously broken within 1 ps (Fig. [2k), yielding OCOC 2 H 4 0~, the majority predicted product 
on Li metal surfaces (Fig. [2b) and a precursor to CO. Here the monovalent anion inter- 
mediate is stabilized by hydrogen bond donation from several AlOH groups and by coordi- 
nation to two surface Li atoms. Since the DFT/PBE method erroneously underestimates 
the 0.51 eV e~ tunneling barrier associated with molecular reorganization (Fig. UK) which 
precedes bond-breaking, it vastly overestimates the overall bond-breaking rate. Indeed, the 
EC decomposition timescale predicted with DFT/PBE is similar similar to the timescale 
predicted in the absence of the ALD layer.— This is in disagreement with our experimental 
measurements which reveals far less solvent decomposition products when an ALD layer is 
present (Sec. [V]). Instead, 0.51 eV should be taken as the overall activation energy in these 
bond-breaking events. With this barrier, the bond-breaking rate should occur in millisec- 
ond, not picosecond, timescales at room temperature. Nevertheless, this PBE-based AIMD 
calculation is valuable because it identifies the most reactive surface site. An EC adsorbed 
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at this site is used in the e~ transfer calculation of the previous section (Fig. H]). Under 
UHV-like conditions, an isolated EC molecule adsorbed at this site exhibits < 0.05 eV adi- 
abatic DFT/PBE C-0 bond-breaking barriers provided a 0.4 V/A electric field is applied 
(Fig. |6]). The qualitative correspondence between adiabatic AIMD/PBE decomposition rate 
and UHV barrierless reaction is the reason this model is adopted for e~ transfer studies in 
Sec. UVBl 

Because of its extreme thinness, optimizing the 7 A -thick LiA102 film coated on to 
Li x C6 has caused 2 Li atoms per surface to migrate outwards (Fig. [2k). These outlying 
Li coordinate to the surface hydroxyl groups, polarizing them. The EC that undergoes 
breakdown (Fig. |6H) is indeed hydrogen bonded to an OH group coordinated to a surface 
Li + . Such Li migration to the surface does not occur in the thicker LiA102 coating. Hence 
the faster adiabatic AIMD /PBE EC decomposition dynamics on the thin LiA102 coating is 
not just a consequence of oxide thickness, but is partly due to active site chemical specificity. 
This anomaly may also be the reason the predicted Vab value does not strongly decrease 
with increasing the oxide thickness from 7 A to 10 A, and may further explain the difference 
in work functions between Li x Cq coated with 7 A and 10 A thick LiA102 films (Table HT|) . 

The 10 A -thick LiA102 coating does not exhibit outward Li atom migration. Here the 
DFT/PBE bond-breaking barriers of adsorbed EC are not readily deconvolved from e~ 
transfer (S.I.). For simplicity, we consider a model with just one 10 A thick L1AIO2 layer 
hydroxylated on both sides (Tabled]), add one excess e~ that now always resides on the EC 
because of the Lia,C 6 e~ sink has been removed, and compute EC~ decomposition energetics 
without applied electric fields. Ce-0 bond-breaking to form COg~ precurors remain barri- 
erless and exothermic. However, the Cc-0 cleavage route to form CO precurors becomes 
endothermic and exhibits a 0.71 eV barrier. This indicates a product channel cross-over 
as the oxide thickness increases and/or the reactivity of the surface site decreases. The 
expected reaction pathyways transition from a mixture of Ce-0 and Cc-0 bond breaking to 
predominantly Ce-0 cleavage (COg - precursor). While the liquid solvent environment is not 
included here, we speculate that this finding may be extrapolated to other coating surfaces, 
including natural SEI films, as the surface sites become less reactive. In the future, we will 
also examine EC decomposition reactions on Li 2 C03 surfaces to see if similar trends per- 
sist on that crystalline material, recently adopted as a theoretical model for organic solvent 
decomposition SEI film, and the decomposition of other solvent/salt molecules.— 
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We have also conducted AIMD simulations of graphitic anodes coated with 5 A thick 
hydroxylated AI2O3 layers (Fig. |2]F). No Li ions reside near the interface region, and no 
solvent decomposition is observed within 7 ps, despite the thinness of the oxide. This 
emphasizes the importance of surface heterogeneity at atomic lengthscales. Replacing all 
surface AlOH groups with AlOLi dramatically increases the decomposition rate; this will be 
discussed in future publications. 

V. EXPERIMENTAL RESULTS 

Figure [7J shows the combined voltammetric and microgravimetric responses of the un- 
coated and alumina coated PLD carbon films as the electrode potential is decreased to a 
value slightly above the threshold for Li + intercalation in the carbon. The uncoated carbon 
electrode (Fig. [Th.) exhibits a continuously increasing current response, with several discrete 
maxima. One maximum reaches a value of 4 //A/cm 2 with a mass increase of 2 /ig/cm 2 at a 
potential of 2 V. The other maximum reaches 11 ///cm 2 at a potential of 1 V. The decompo- 
sition of the electrolyte and deposition of byproducts at 2 V is catalyzed by the Cu substrate, 
as evidenced by the similar current and mass changes on a control Cu electrode (Fig. [7b), 
and demonstrate that the carbon films possess porosity and allow electrolyte penetration. 
As seen in the limiting current and mass profiles of Fig. [7b, electrolyte decomposition results 
in Cu passivation beyond 2 V, arguing that the majority of the current and mass changes 
measured above 2 V for the porous carbon films (Fig. [TJi) are due to electrolyte decom- 
position on Cu; only the signal below 1 V is associated with solid electrolyte interphase 
formation (SEI) on the carbon surface. The porous and therefore higher area carbon surface 
exhibits a continuous increase in both current and mass uptake as the potential is further 
reduced from 2 to 0.2 V and the onset of Li + intercalation is approached. With an approach 
to 0.2 V, the rate of current change increases substantially over the rate of mass change, 
signaling a point where Li + intercalation has initiated, where the lighter mass Li (compared 
to an fragment of ethylene carbonate or diethyl carbonate) accounts for a growing fraction 
of the measured current. The possibility exists that current increase could also be related to 
solvent reduction without mass addition to the surface (soluble byproduct formation), but 
note that the mass decrease upon reversal of the potential sweep clearly argues for the onset 
of Li + ion intercalation into the carbon. We note that the scan rate of 1 mV/s is sufficiently 
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fast to produce only modest extraction of Li during this reverse partial half cycle. 

The alumina coating acts as a kinetic barrier to prevent electron transfer to the organic 
carbonate molecules of the electrolyte. Figures [7b, d show the response of a 0.55 and 1.1 nm 
thick coated carbon films to the onset of electrolyte reductive decomposition. Comparison 
of the uncoated (Cu subtracted) and coated carbon films shows that a higher overpoten- 
tial is required to drive solvent decomposition and a lower quantity of mass addition takes 
place with the alumina coating present. A Cu current and mass uptake response is elimi- 
nated for these coated electrodes because the alumina nucleating agent and film precursors 
fully penetrate the porous carbon, conformally coating both the carbon network and the 
underlying exposed regions of the Cu substrate. The onset for significant current density 
and mass increase occurs at approximately 1.2 V and 0.8 V for the 0.55 nm and 1.1 nm 
alumina coatings, compared to 1.5 V for the uncoated carbon. Mass increases measured 
at 0.8 V are 6, 1.3 and 0.5 /zg/cm 2 for the uncoated (Cu subtracted), 0.55 nm and 1.1 nm 
alumina sample, respectively. The greater overpotential and reduced mass uptake of the 
1.1 nm coating relative to the thinner 0.55 nm coating argue that the thicker film provides a 
more effective kinetic barrier for reducing the extent of both reductive solvent decomposition 
and byproduct deposition on the electrode. The thicker alumina film would be expected to 
present a lower electron tunneling rate resulting in a slower rate of solvent decomposition 
and retarded SEI formation. The fact that mass addition is observed in the presence of these 
alumina coatings is a clear indicator that alumina serves to retard and limit the extent of 
but does not prevent electrolyte reduction and resulting byproduct film formation. 



VI. CONCLUSIONS 



In this work, we compare EC decomposition on Li metal and on models of oxide-coated 
electrodes. The latter mimics recent experimental work using ALD technique to passivate 
anodes. This ALD strategy carries significant technological promise,- - — and it also provides 
an ideal robust platform for theoretical and experimental study of passivating mechanisms. 
These two systems represent two electron transfer regimes. 

On pristine Li (100) surfaces, liquid EC and even isolated adsorbed EC molecules are 
predicted to undergo decomposition in picosecond time scales. CO is the dominant product 
from EC, possibly because of favorable kinetic prefactors, even though both the CO and 



25 



COg - reaction pathways are almost barrierless and the COg - product is more thermody- 
namically stable. EC molecules and the electrode are in close contact and strongly coupled. 
Adiabatic DFT/PBE and AIMD/PBE simulations should be accurate in this regime. 

In contrast, electron transfer through an oxide layer should be slow compared to nuclear 
motion. We find evidence that tunneling through even a 7 A thick oxide layer belongs 
to the non-adiabatic regime. Applying constrained DFT (cDFT) calculations, such thin 
coatings are found to slow down e~ transfer because the solvent reorganization energy A now 
figures prominently in electron tunneling through the oxide. A, largely neglected in previous 
studies of electrolyte decomposition in batteries, is estimated to be ~ 2 eV for adsorbed 
EC molecules in ultra-high vacuum-like conditions. This translates into a ~ 0.5 eV electron 
tunneling barrier within the harmonic approximation when the e~ transfer free energy change 
is small. 

cDFT calculations show that the 7 A - and 10 A -thick LiA102 coated Li^Ce exhibit elec- 
tron transfer rates of ~ 10 5 /s at the Li + /Li(s) applied voltage. The predicted e~ transfer 
rate is not free of ambiguities and assumptions, and is of order-of-magnitude utility; fur- 
ther fundamental research is needed for a more rigorous treatment. Despite this caveat, 
this work respresents the first first-principles estimate of the e~ tunneling rate between an 
electrode and an EC molecule across an insulating oxide layer. Such predictions are critical 
for understanding ALD-hindered SEI growth in lithium ion batteries. 

The overall electron transfer rate (Eq. [T]or Eq. [5]) also depends on the offset AG Q between 
e~ donor and acceptor species. AG Q in turn depends on the applied voltage. AIMD estimates 
of AG Q in an explicit liquid solvent environment is currently lacking, and we have relied on 
dielectric continuum treatments of the liquid environment. Nevertheless, our analysis yields 
useful insights. With any reasonable estimate of AG Q , the electron transfer rate to EC at 
the surface is predicted to be faster than 1/s, and solvent breakdown on the ALD oxide 
is expected. This is confirmed by our gravimetric measurements on ALD-coated anodes, 
although the amount of solvent decomposition product is significantly less than that on 
uncoated graphite electrodes. 

In the case of oxide-coated electrodes, AIMD/PBE and DFT/PBE calculations with- 
out electronic constraints vastly underestimate the electron transfer barrier. The reason 
is most likely the self-interaction error, which unphysically favors a split electron partially 
localized on the EC and partially delocalized on the electrode. This defect exists in many 
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DFT functional and has been known to yield errors in when a molecular is split into two 
fragments. 32 As a result, direct AIMD/PBE simulations overestimate EC decomposition 
rates at oxide-coated electrode surfaces by many orders of magnitude. However, AIMD/PBE 
and DFT/PBE calculations still provide a wealth of information about structure and rela- 
tive energetics, and they form the basis of Marcus theory considerations and non-adiabatic 
electron transfer studies which are key aspects of this work. 

Taking advantage of the qualitative correspondence between AIMD liquid state reaction 
rates and ultra-high vacuum-like DFT calculations of barrier heights at T=0 K in an electric 
field, we have applied calculations in UHV-like settings to suggest that the dominant product 
may shift from a mixture of CO and CO 2- to mainly CO 2- as the binding of ionic decom- 
position products becomes less favorable (e.g., on thicker oxide coatings). This prediction 
may be transferrable to natural SEI films arising entirely from electrolyte decomposition. 
Atomic-scale surface heterogeneity is found to affect EC decomposition, with Li + ions at the 
surface playing a facilitating or "catalytic" role. Our work paves the way for novel future 
experimental studies in UHV settings. 
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(b) 




FIG. 1: (a) Schematic of electron transfer between isolated orbitals. The red and blue represent 
the diabatic potential energy surfaces of e~ donor and acceptor as a function of the generalized 
polarization (P) degree of freedom. The green segments represent adiabatic processes with the 
non-crossing surfaces split by 2Vab- AG is the reaction free energy and A is the reorganization 
energy. The green roman numbers denote (i) flat EC; (ii) flat EC~; (iii) bent EC - ; (iv) bent EC. 
(b) Non-adiabatic e~ transfer form a metallic electrode. The thick upper red line represents the 
Fermi level, and is the primary donor orbital within cDFT calculations. e~ can also transfer from 
the continuum of electrode donor states below the Fermi level, depicted as thin red lines, to the 
acceptor orbital, with however increased non-adiabatic barriers (crossing points between blue and 
red curves), (c) & (d) Flat and bent EC molecules, respectively. Red, grey, and white refer to O, 
C, and H atoms. 



34 




FIG. 2: (a) Model electrode system with a narrow LiCq strip decorated with C=0 edges and coated 
with a 10 A thick LiA102 layers terminated by hydroxyl groups, (b) A snapshot 15 ps into an AIMD 
trajectory of 32 EC molecules confined between Li metal (100) slabs conducted at T=350 K. Li and 
decomposed EC are depicted as ball-and-stick models, intact EC as wireframes, (c) An isolated, 
decomposed EC on Li surface at T=350 K, after a 7 ps AIMD simulation, (d) No decomposition 
on 10 A thick LiA102-coating after 7 ps. (e) One EC decomposed on 7 A thick LiAlOvcoated 
surface, T=450 K, after 7 ps. (f) No decomposition on 5 A thick, hydroxylated A^Os-coating 
after 7 ps. Yellow, grey, red, blue, and white depict Al, C, O, Li, H atoms, respectively. 
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FIG. 3: (a) Static, nudged elastic band (NEB) calculations of energy barriers associated with 2 
modes of EC breakdown, producing CO or CO3 - , at T=0 K. Points B, C, & D correspond to panels 
(b)-(d). (b) Intact EC on Li (100). (c) EC partially decomposed into OCOC2H4O, precursor to 
CO and OC2H4O 2- , on Li (100). (d) COg~ and C2H4 products on Li metal. The color key is the 
same as in Fig. [2j 
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FIG. 4: (a)&(b): Local electronic densities-of-state decomposed along the x-axis (perpendicular 
to interface) for EC adsorbed on thin LiA102-coated LiCg. Panels (a) and (b) correspond to the 
configurations depicted in panels (c) (flat EC geometry) &(d) (bent geometry), respectively. The 
red patches depict integrated up- and down-spin densities exceeding 0.01 |e| for each planewave 
wavefunction collapsed on an atom centered at x. Panel (b) shows that the bent geometry dras- 
tically changes the HOMO and LUMO levels, with an excess electron now residing on EC below 
the Fermi level (-Ef)- The conduction band of the LiA102 region is located above 4 eV. In panels 
(c) & (d), the EC molecule and Li coordinated to the EC are depicted as spheres while surface 
hydroxyl groups donating hydrogen bonds to the EC are stick figures. Other Li are omitted and 
all other oxide-coating and graphite atoms appear as wireframes. EC configurations on the thick 
LiA102 coating are qualitatively similar (not shown). The color scheme is as in Fig. [2j 
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FIG. 5: (a)&(b) Adiabatic (red crosses) and non-adiabatic (blue crosses) energy profiles along 
the reaction coordinate between the flat (Fig. [3b, image "0") and bent (Fig. 0H, image "5") EC 
geometries when applying a 0.4 V/A applied electric field to the 7 A and 10 A-thick L1AI2O 
layers. Adiabatic energies are computed along the NEB-generated chain with unconstrained DFT. 
Non-adiabatic reorganization energies (A) derive from cDFT. The dashed curves are parabolic fits. 
(c)&(d) Highest occupied orbital of the system (green), and the cDFT-computed e~-accepting 
EC orbital (red) adsorbed on the 7 A and 10 A-thick L1AI2O layers, integrated over the lateral 
dimensions. 
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FIG. 6: (a) Adiabatic DFT/PBE energy profiles associated with EC decomposition on the thin 
LiAlC-2 coating at T=0 K. Red and blue refer to adsorbed OCOC 2 H 4 2 " and C2H4OCO2" inter- 
mediates, which are precursors to CO and CO 2- products, respectively. The electric field strength 
is 0.4 V/A. The dashed lines are guides to the eye. Red/blue triangles depict the barrier associated 
with the first transfer of an electron on to the EC molecule, which is detected between images. 
They indicate that the DFT/PBE treatment erroneously neglects the electron tunneling barrier. 
Black arrows denote first detection of C-0 bond breaking in an image. Point B correspond to the 
intact EC in Fig. 3b. (b) Same as (a) but for a single 10 A thick layer in the absence of Li^Cg or 
electric field, (c)-(d) OCOC2H4O 2- , and C2H40C02~ on thin LiA10"2 surfaces, corresponding to 
points C & D in panel (a). The atom representation is as in Fig. [U 
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FIG. 7: Current and mass change response with 1 mVs cathodic polarization of an electrode in 
1 M L1PF6, 1:1 vol. EC:DEC. (a) 50 nm C film, (b) Cu substrate without a C film, (c) 0.55 nm 
thick ALD AI2O3 on a 50 nm C film, and (d) 1.1 nm thick ALD AI2O3 on a 50 nm C film. Current 
and mass are normalized to the geometric area. 
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